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We present an in-depth analysis of the spectral and temporal behavior of a long almost uninter- 
rupted INTEGRAL observation of Vela X-1 in Nov/Dec 2003. In addition to an already high 
activity level, Vela X-1 exhibited several very intense flares with a maximum intensity of more 
than 5 Crab in the 20^10 keV band. Furthermore Vela X-1 exhibited several off states where 
the source became undetectable with ISGRI. We interpret flares and off states as being due to 
the strongly structured wind of the optical companion: when Vela X-1 encounters a cavity in the 
wind with strongly reduced density, the flux drops, thus potentially triggering the onset of the 
propeller effect which inhibits further accretion, thus giving rise to the off states. The required 
drop in density to trigger the propeller effect in Vela X-1 is of the same order as predicted by 
theoretical papers for the densities in the OB star winds. The same structured wind can give rise 
to the giant flares when Vela X-1 encounters a dense blob in the wind. Further temporal analysis 
reveals that a short lived QPO with a period of 6800 sec is present. The part of the light curve 
during which the QPO is present is very close to the off states and just following a high intensity 
state, thus showing that all these phenomena are related. 
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1. Introduction 

Vela X-l is a high mass X-ray binary (HMXB) consisting of the super giant HD 77581 and 
a massive (1.9 M Q [1]) neutron star in a 8.964 day orbit [2]. The optical companion has a mass 
of ~23M and a radius of ~30/? Q [2]. The neutron star is deeply embedded in the dense stellar 
wind of the donor HD 77581 (M (c = 4x 1CT 6 M yr" 1 ) [3]. X-ray lines indicate that this wind is 
inhomogeneous with many dense clumps [4] embedded in a far thinner, highly ionized component 
[5]. 

The neutron star has a long spin period of ~283 s [6]. The evolution of the spin period is 
best described by a random walk as expected for a wind-accreting system [7]. Although the source 
exhibits strong pulse-to-pulse variations, a pulse-profile folded over several pulse periods shows 
remarkable stability [8], even over decades [9]. At energies below 5 keV, the pulse-profile consists 
of a complex five-peaked structure, which transforms at energies above 20 keV into a simple double- 
peaked pulse-profile [8] where the two peaks are thought to be due to the two accreting magnetic 
poles of the neutron star. 

With an X-ray luminosity of ~4 x 10 36 erg s _1 , Vela X-l is a typical high mass X-ray binary. 
Previous observations have shown that the source is strongly variable with reductions to less than 
10% of its normal value [10, 11, 12, 13], while periods of increased activity have also been observed 
during which the flux increases within an hour to a multiple of the previous value, reaching peak 
flux levels close to lCrab [11, 14, 15]. In this respect, Vela X-l is similar to sources such as 
4U 1700—377 and 4U 1907+09, for which low luminosity states and flares have also been observed, 
as is rather typical for wind-accreting systems [16, 17, 18, 19]. Although Vela X-l is a well studied 
object, only observations by INTEGRAL revealed that the flares in Vela X-l can be brighter than 
previously anticipated [10, 20, 21]. 

2. Data 

INTEGRAL observed the Vela region continuously for five consecutive INTEGRAL revo- 
lutions from revolution 137 (JD 2452970.86) until revolution 141 (JD 2452985.44) resulting in 
approximately 1 Msec of data (see Fig. [I]). 

We used all available science windows (ScWs) to be able to derive a contiguous light curve 
with as few interruptions as possible. Since Vela X- 1 is a bright source, the OS A can detect the 
source and determine its flux level accurately even when the source is at high off-axis angles. 

Apart from Vela X-l, also 4U 0836-429, H0918-5459, the Vela Pulsar, and two sources first 
reported by INTEGRAL [22, 23] are detected. Vela X-l is widely separated from the other sources 
such that contamination of the spectrum of Vela X-l is of no concern. Data from JEM-X and SPI 
have not been used in this analysis due to the far smaller field-of-view of JEM-X and since Vela X-l 
is off-center in the observed field (Vela X-l was only within the fully coded field-of-view of JEM-X 
for less than ten out of the ~550 individual pointings). 

3. Data analysis 

Vela X-l was found in a highly variable state during the observation. While Vela X-l is known 
to be a variable source [11, 24], the behavior found in this observation [10, 20] is indeed extreme. 
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Figure 1: a) ISGRI ScW by ScW light curve in the 20^0 keV energy band and b) hardness ratio. The labels 
indicate the revolution number (from [10]). 



Most importantly, on 2003 November 28 (JD 2452971.67), INTEGRAL observed an extremely 
bright flare (flare 1; see Fig. [j]). During the flare, the 20^-0 keV count rate increased from a ScW 
averaged pre-flare value of ~55countss~ 1 (~300mCrab, or 1.6 x 10~ 9 ergcm~ 2 s _1 ) by a factor 
of more than seven to 405 counts s _1 (2.3 Crab) within only 90 minutes. 

After the peak, the flare decayed quickly to an intensity level of <1 Crab and within ~11 h to 
a ScW averaged count rate of ~35 counts s _1 (200mCrab), somewhat lower than before the onset 
of the flare (see Fig. II]). 
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Figure 2: Close up of the light curve of giant flare 1 with a time resolution of 20 s. The peak is reached at 
MJD 52971.18 with 923 countss" 1 (corresponding to 5.2 Crab) in the 20 to 40keV band (from [10]). 



Table 1: Overview of the observed flares. See Fig. 1 for the numbering of the flares. The time is the onset 
of the flare. To obtain the peak fluxes F p eak> a light curve with a time resolution of 20 s was used. r r j se is the 
time from the onset of the flare to the peak, while T tota \ is the duration of the flare. 



Flare 


Time 
[MJD] 


Duration 
[s] 


^peak 

[Crab] 


T- 

1 rise 
^total 


Remarks 


1 


52971.15 


11200 


5.2 


0.15 


giant flare, spectral softening 


2 


52975.34 


5 200 


2.6 


0.83 


no spectral change 


3 


52976.50 


1800 


5.3 


0.28 


giant flare, very short 


4 


52977.15 


12 900 


1.9 


0.13 


spectral softening 


5 


52980.31 


31400 


3.9 


0.63 


high intensity state, no spectral change 



In the following 13 days, three more flares (flares 2 to 4, see Table [I] were observed. All 
three flares were shorter and less intense than flare 1 on a science window averaged basis, but still 
reached ScW averaged intensities close to 1 Crab. 

On 2003 December 7 (JD 245981.10), another intense flare was observed (designated Flare 5, 
see Fig. |]). Unlike flare 1, during which the brightness of the source increased rapidly, it took ~8 h 
for flare 5 to reach its ScW averaged maximum 20-40 keV flux of ~ 1.2 Crab. The decay lasted 
~5h until Vela X-l reached its pre-flare count rate of ~35 counts s _1 (200mCrab in 20-40 keV). 
Although quite bright, flare 5 is therefore significantly less intense than giant flare 1, and also far 
longer, i.e. it is a high intensity state. 

The analysis of a light curve with a 20 s time resolution showed that the source reached a 
peak count rate of ~920 counts s _1 (5.2 Crab) in flare 1 and ~930 counts s _1 (5.3 Crab) in flare 3. 
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Figure 3: Closeup of the light curve (with a time resolution of 283 s to remove the pulsations) where the 
temporary QPO is present. Note that during the trough between 2 h and 3 h, and especially following the 
quasi-periodic modulation, the count rate decreased several times to zero for a short time (from [10]). 

Flare 3 on December 3 was therefore also a giant flare. However, flare 3 was significantly shorter: 
the entire flare lasted less than 2000 s, but it was as bright as flare 1 (see Table |]). 

Extending the analysis to the non-flaring parts of the light curve, we detected a quasi-periodic 
oscillation (QPO), similar to other accreting X-ray pulsars. The short-lived QPO with a period of 
~6820 sec appears to be quite regular and inconsistent with pure stochastic behavior (see Fig. ||). 
Subsequent period searches on the corresponding data subset clearly detect the period. We note that 
the quasi-periodic modulation shown in Fig. || is far stronger and inconsistent with the NOMEX 
effect, which can cause intensity variations from pointing to pointing, but not within a given point- 
ing. 

Furthermore, we observed several off states, during which no significant residual flux was 
detectable by ISGRI (Fig. Q). The onset of these off states occured very suddenly. The luminosity 
of the source simply drops below the detection limit of ISGRI. At the end of the off states, Vela X-l 
switches instantly on again and immediately resumes its normal intensity level. All off states 
occured within 12h from MJD 52981.0 and MJD 52981.5. 



4. Spectral evolution during the flares 

Although the source was extremely bright during the flares, meaningful spectral fits could 
not be obtained, as the exposure time was to short. We therefore analyze the hardness ratio (see 
Fig. [j]b), which is defined as 

H-S 

HR = (4.1) 

H + S 
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Figure 4: Close-up with a time resolution of 20 sec on off states 3 and 4 of Vela X-l during which the source 
becomes undetectable by ISGRI and then turns on again within one hour (from [10]). 



where H is the count rate in the hard band (40-60 keV) and S the count rate in the soft band 
(20-30 keV). While the hardness ratio remained constant throughout most of the observation at 
~ —0.735, the hardness ratio significantly changed with the onset of flare 1: it dropped to —0.82 
and during the flare to —0.85 (see Fig. |l|b). The same for flare 4: the hardness ratio dropped from 
—0.72 to —0.84, the same level as in giant flare 1, although flare 4 was far shorter and reached only 
a third of the peak flux of flare 1. During flares 2 and 5, however, the hardness ratio did not change. 

We then used "Hardness intensity diagrams" (HIDs) to study the spectral evolution of the 
source (see Fig. ||). Most of the data points are centered around the average values of intensity and 
hardness ratio. The only exception are the data points from the flares, which are above the general 
cluster of data points. Due to the spectral softening, the data points of flare 1 and flare 4 are shifted. 
The softening, however, did not evolve during the flares, but the flares are softer than the average 
spectrum from the beginning until the end. 



5. Discussion 
5.1 The flares 

Vela X-l has always been known to be highly variable with time and to show intensity varia- 
tions of up to a multiple or a fraction of the original intensity on all time scales. Although Vela X-l 
has exhibited extensive flaring activity in the past however, giant flares (as flares 1 and 3) had not 
been seen before. 

The analysis of the hardness ratio shows that there seem to be two types of flares: the first type 
(flares 1 and 4) shows dramatic increases in the count rate, the onset of the flare is very sudden, 
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Figure 5: Hardness Intensity Diagram of Vela X-l; data from the eclipses have been excluded. The data 
points from the flares are indicated by individual symbols. The datapoints from the eclipses have been 
removed (from [10]). 



and the spectrum softens during the flare. The second type is similar to a high intensity state: these 
flares are longer and the spectrum does not change. 

The mechanism behind the flaring activity, however, is not fully understood. It has been shown 
that a temporary disk may form in Vela X-l [25]. The disk collapses and the material is accreted 
onto the neutron star giving rise to a short flare. These predicted flares would last from 15 to 60 
minutes, similar to the short flaring activity of Vela X-l. Furthermore, wind accretion is a highly 
instable process by itself: the accretion wake trailing the neutron star contains filaments which also 
produce flares when being accreted [26]. The shock trailing the neutron star oscillates creating the 
"flip-flop instability" which then produces inflows that repeatedly change their direction [27, 28]. 
The timescale of 45 min matches some of the observed behavior very well, but fails to explain the 
long flares. 

As the local density in a shocked wind varies by a factor of 100 [29], which can explain the 
flaring X-ray luminosity [4], dense clumps trapped in an otherwise thin and more homogeneous 
wind might be responsible for long flares [30], when the clumps are being accreted. Such a clump 
can feed the neutron star with a significantly higher M than usual over several hours. In summary, 
we conclude that the observed long flares are due to a strongly structured OB star wind, while when 
Vela X- 1 is less active, the OB star wind is less structured. 

5.2 The off states 

In a similar way to the flaring activity, the off states [13, 11, 10] where the source is below 
the detection limit are remarkable. After the off state observed by RXTE in 1996 [11], the source 
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resumed its normal, pulsating behavior without any transition phase. The offstates reported here 
also occured without a transition phase [10]. 

The reasons for these off states and the sudden reappearance of pulsations are not understood. 
So far, several ideas to explain these phenomena have been proposed, however, none can fully 
explain the observed off states, since all of these ideas require a significantly longer timescale, e.g. 
clumps in the stellar wind [31] would have to have an unrealistic high optical depth to completely 
block hard X-rays and can not pass the line of sight in a few seconds to explain the sudden turn 
on/off behavior of the source. Therefore other mechanisms must be considered. 

The wind of OB super giants is inhomogeneous and clumpy [32, 26]. The density of the 
stellar wind can vary by several orders of magnitude [33]. Not only clumps, but also holes, i.e. 
regions of strongly reduced density are present in the wind: there, the density is lower than the 
average density of the wind by a factor of 10 3 [33]. If the neutron star enters these holes, M would 
then also decrease by a factor of ~ 10 3 and the X-ray luminosity would be reduced accordingly. 
Furthermore, the density fluctuations predicted by these models occur suddenly [see Fig. 1 in 33] 
similar to the onset of the off states (see Fig. |]). If M drops due to these density variations in the 
wind, the Alfven radius will increase due to the reduced ram pressure of the infalling gas. Once it is 
larger than the co-rotation radius, accretion onto the neutron star is inhibited, i.e. the X-ray source 
basically switches off. This scenario is commonly known as the propeller effect [34]. Since the 
propeller effect depends on the amount of infalling material, the Alfven radius is not constant. This 
effect was observed in GX 1+4 [35]: in very low luminosity states, no pulsations were observable, 
while the source was strongly pulsating in high luminosity states. 

Since the strength of the magnetic field of Vela X-l is known from the observation of the 
cyclotron lines [36], the critical flux limit for Vela X-l for the onset of the propeller effect can 
be obtainted (after [35]): ^x,Propeiier,VeiaX-i w 1.1 x 10~ 12 erg cm _2 s _1 . Compared with the typical 
bolometric flux of several times 10~ 9 erg cm s , this critical flux is lower by about three orders 
of magnitude. This flux limit matches very closely the predicted density variations in the stellar 
wind of 10 3 5 [32]. 

We therefore conclude that off states could be caused by a sudden drop in M that allows 
Vela X-l to enter the propeller regime. Intensity dips, however, are longer, show a smooth transi- 
tion, and exhibit photoelectric absorption of more than 10 24 cm~ 2 . These dips are readily explained 
by a dense blob in the wind passing through the line of sight. 

5.3 Connection with SFXTs 

The similarity between the flares and off states in Vela X-l and the behavior of Supergiant 
Fast X-ray Transients (SFXTs) [37] is striking. SFXTs are high mass X-ray binaries that show very 
brief outbursts on timescales of hours or even only tens of minutes, and then remain undetectable at 
higher energies for months between outbursts [38]. SFXTs should be rather bright persistent objects 
[39] since the neutron star is deeply embedded in the dense stellar wind of the optical companion, 
however, the accretion is inhibited by the propeller effect. It has therefore been proposed that 
SFXTs harbor a magnetar [40] and that the extremely strong magnetic field (B > 10 15 G) effectively 
inhibts accretion unless the ram pressure of the infalling gas is high enough such that accretion 
becomes possible for a short time. The giant flares and off states of Vela X- 1 are therefore similar 
to these outbursts: in both cases the accretion of a dense blob of material causes the outburst or 
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flare, while reduced material infall causes the source to switch off. In summary, Vela X-l and 
SFXTs are rather similar objects, however, SFXTs are usually in the off state, while Vela X-l is 
usually in a normal accretion mode. 
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